Persoonia 41, 2018: 56-70 
www.ingentaconnect.com/content/nhn/pimj 


ISSN (Online) 1878-9080 | 
https://doi.org/10.3767/persoonia.2018.41.04 


RESEARCH ARTICLE 


Malassezia vespertilionis sp. nov.: 
a new cold-tolerant species of yeast isolated from bats 


J.M. Lorch’, J.M. Palmer?, K.J. Vandervvolf 15, K.Z. Schmidt’, M.L. Verant'?, 
T.J. Weller^, D.S. Blehert' 


Key words Abstract Malassezia is a genus of medically-important, lipid-dependent yeasts that live on the skin of warm- 
blooded animals. The 17 described species have been documented primarily on humans and domestic animals, 
but few studies have examined Malassezia species associated with more diverse host groups such as wildlife. 
While investigating the skin mycobiota of healthy bats, we isolated a Malassezia sp. that exhibited only up to 92 % 
identity with other known species in the genus for the portion of the DNA sequence of the internal transcribed spacer 
region that could be confidently aligned. The Malassezia sp. was cultured from the skin of nine species of bats 
in the subfamily Myotinae; isolates originated from bats sampled in both the eastern and western United States. 
Physiological features and molecular characterisation at seven additional loci (D1/D2 region of 26S rDNA, 18S 
rDNA, chitin synthase, second largest subunit of RNA polymerase II, B-tubulin, translation elongation factor EF-1a, 
and minichromosome maintenance complex component 7) indicated that all of the bat Malassezia isolates likely 
represented a single species distinct from other named taxa. Of particular note was the ability of the Malassezia 
sp. to grow over a broad range of temperatures (7—40 °C), with optimal growth occurring at 24 °C. These thermal 
growth ranges, unique among the described Malassezia, may be an adaptation by the fungus to survive on bats 
during both the host's hibernation and active seasons. The combination of genetic and physiological differences 
provided compelling evidence that this lipid-dependent yeast represents a novel species described herein as 
Malassezia vespertilionis sp. nov. Whole genome sequencing placed the new species as a basal member of the 
clade containing the species M. furfur, M. japonica, M. obtusa, and M. yamatoensis. The genetic and physiological 
uniqueness of Malassezia vespertilionis among its closest relatives may make it important in future research to 
better understand the evolution, life history, and pathogenicity of the Malassezia yeasts. 
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INTRODUCTION has even been hypothesised that Malassezia could play a role 
in promoting certain forms of skin cancer (Gaitanis et al. 2011). 
Members of the genus Malassezia are lipid-dependent fungi 
specialised to live on the skin of humans and other euthermic 
animals. Malassezia is the sole genus in the class Malassezio- 
mycetes, and the 17 described species appear to be part of the 
natural skin mycobiome of animals (Wang et al. 2014). At least 
five species have been regularly associated with dermatitis or 
other types of skin disorders in humans (reviewed by Gaitanis 
et al. 2012), and M. pachydermatis is described as the cause 
of otitis externa in domestic dogs (Gustafson 1955, Bond et 
al. 2004). However, these potentially pathogenic species of 
Malassezia are also found on areas of normal skin of afflicted 
patients and on asymptomatic individuals, making it unclear 
what role the fungi play in skin disease (reviewed by Gaitanis 
etal. 2012). Malassezia furfur, M. pachydermatis, and M. sym- 
podialis have also occasionally been implicated as causes of 
sepsis in infants and immunocompromised patients (reviewed 
by Gaitanis et al. 2012, Aguirre et al. 2015, Patron 2016), and it 


In addition to their medical ambiguity, relatively little is known 
about the diversity and ecology of Malassezia. Members of 
the genus were traditionally identified on the basis of pheno- 
typic traits, and prior to 1996 there were only three recognised 
species of Malassezia. More recent application of molecular 
techniques to assist with species characterisation has facilitated 
the ability to distinguish cryptic species and has increased the 
known diversity of the genus (e.g., Sugita et al. 2002, Hirai et 
al. 2004, Cabafies et al. 2007, 2016, Honnavar et al. 2016). 
Sugita et al. (2010) present a list of the animal hosts from 
which various species of Malassezia have been recovered. 
However, many reports used to generate that list lacked 
molecular data to support the identification of the Malassezia 
species that were isolated, and potentially novel taxa may have 
been overlooked. Indeed, diversity of the genus is likely much 
higher than currently documented (Amend 2014, Cabafies 
2014). Ten of the 17 Malassezia species are most closely as- 
sociated with humans, and were discovered through culture- 
based surveys of diseased skin; the remaining seven species 
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(M. brasiliensis, M. caprae, M. cuniculi, M. equina, M. nana, 
M. pachydermatis, and M. psittaci) have been isolated primarily 
from animals (reviewed by Sugita et al. 2010, Cabafies et al. 
2016), with M. pachydermatis reported as a zoonotic pathogen 
(Chang et al. 1998). Some of the zoophilic members of the 
group appear to have a broad host range, while others are more 
host-specific (reviewed by Guého-Kellerman et al. 2010, Sugita 
et al. 2010). Little work has been done with broader taxonomic 
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host groups, and the relatively small number of described spe- 
cies of Malassezia is likely the result of sampling bias, which 
is skewed toward humans and domestic animals. Given the 
host specificity of some species of Malassezia, many more 
taxa may be discovered when a broader range of host species 
(especially wildlife) are sampled. Such undiscovered species 
of Malassezia could be important in further elucidating the 
taxonomy, evolution, ecology, and pathogenicity of this group 
of medically important fungi. 


While investigating the mycobiota on the skin of bats, we de- 
tected a putative Malassezia sp. that was genetically distinct 
from other known members of the genus. Here we describe the 
isolation, occurrence, and characterisation of this novel species. 


MATERIALS AND METHODS 


Isolation of Malassezia from bats 


Samples were collected in 2014, 2016, and 2017 underthe U.S. 
Geological Survey - National Wildlife Health Center (NWHC) 
Animal Care and Use Committee Protocols #EP140212 and 
#EP081124-A2, with all necessary permits and permissions 
for the sites and species sampled. Hibernating bats (Fig. 1) 
were captured by hand and active bats were captured in mist 
nets. Gloves were changed between animals to prevent cross- 
contamination. The animals were sampled non-lethally using 
sterile Pur-Wraps® polyester-tipped swabs (Puritan Medical 
Products Company LLC, Guilford, Maine, USA) pre-moistened 
with 150 uL of sterile nuclease-free water. Swabs were gently 
rolled back-and-forth three times across the skin of the fore- 
arm and wing membrane between the elbow and wrist joints. 
Samples were then placed in sterile microcentrifuge tubes, 
stored chilled for up to 48 h, and shipped on ice to the NWHC. 
A total of 264 samples were obtained from thirteen sites in 
seven states (one site in Alabama, USA; one site in California, 
USA; one site in Kentucky, USA; one site in Missouri, USA; 
one site in Pennsylvania, USA; two sites in New York, USA; 
and six sites in Wisconsin, USA), representing ten bat species 
(Lasionycteris noctivagans, Myotis californicus, Myotis grises- 
cens, Myotis leibii, Myotis lucifugus, Myotis septentrionalis, 
Myotis sodalis, Myotis thysanodes, Myotis yumanensis, and 
Perimyotis subflavus). A list of all individual bats sampled for 
the project is provided in Table 1. 
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Upon arrival at the laboratory, swabs were streaked onto Leem- 
ing and Notman agar (LNA; 10 g bacteriological peptone, 0.1 g 
yeast extract, 5 g glucose, 8 g desiccated ox bile, 1 mL glycerol, 
0.5 g glycerol monostearate, 0.5 g Tween 60, 10 mL whole fat 
cow's milk, 0.5 g chloramphenicol, 0.5 g cycloheximide, 15 g 
agarose per litre, pH 6.0; modified slightly from Leeming & Not- 
man (1987)) and incubated at 7 °C. Plates were checked weekly 
for a total of 12 wk, and any colonies resembling Malassezia 
were transferred to fresh LNA. Isolates were identified by se- 
quencing the ITS as described by Lorch et al. (2015). 


Whole genome sequence analysis 


Isolate CBS 15041 (NWHC 44797-103; UAMH 11924) was select- 
ed for whole genome sequencing to further resolve the taxono- 
my ofthe bat-associated Malassezia. Nucleic acid was obtained 
using a phenol-chloroform extraction. Library preparation and 
next-generation sequencing was performed by the University 
of Wisconsin Biotechnology Center DNA Sequencing Facility 
using the genomic Nextera XT DNA Library Prep Kit (Illumina 
Inc., San Diego, CA) and the Illumina MiSeq Next Genera- 
tion Sequencer platform. Sequence data was processed and 
assembled using JAAWS (https://github.com/nextgenusfs/ 
jaaws). Briefly, the paired-end 250-bp MiSeq sequence reads 
(2x250) were processed with trimmomatic v. 0.36 (Bolger et 
al. 2014) to remove adapter sequences and phiX spike-in was 
removed using bowtie2 v. 2.3.2 (Langmead & Salzberg 2012) 
alignment to the phiX genome (NC. 001422). The data were 
then assembled into scaffolds with Spades v. 3.9.0 (Bankevich 
et al. 2012). The subsequent assembly was cleaned using 
Blobtools v. 0.9.19 (Laetsch & Blaxter 2017) and filtered for 
unexpected coverage, mitochondrial DNA, contamination, and 
scaffolds less than 1 kb in length. Finally, the cleaned assembly 
was error corrected using five iterations of Pilon v. 1.22 (Walker 
etal. 2014). The genome of the bat-associated Malassezia was 
annotated with funannotate v. 0.7.0 while the 28 genomes of 
Malassezia species previously sequenced (Wu et al. 2015; 
Table 2) were annotated using funannotate v. 0.5.3 (https:// 
github.com/nextgenusfs/funannotate). Conserved orthologues 
were identified using BUSCO2 (Simao et al. 2015) basidi- 
omycete database using the busco wrapper script in Phyloma 
(https://github.com/nextgenusfs/phyloma). The concatenated 
protein sequences of 254 conserved BUSCO2 orthologues 
were aligned using MAFFT v. 7.305b (Katoh & Standley 2013) 


Fig. 1 Hibernating bats, such as these Myotis sp., were sampled for this study by swabbing wing skin. 
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Table1 List of individual bats sampled for this study. 


Individual identifier 


Location 


Host species 


Sampling date 


Malassezia vespertilionis isolated 


24716-001 
24716-002 
24716-003 
24716-004 
24716-005 
24716-006 
24716-007 
24716-008 
24716-009 
24716-010 
24716-011 
24716-012 
24716-013 
24716-014 
24716-015 
24716-016 
24716-017 
24716-018 
24716-019 
24716-020 
24738-002 
24738-005 
24738-006 
24738-008 
24738-011 
24738-013 
24738-022 
24738-025 
44767-001 
44767-002 
44767-003 
44767-004 
44767-005 
44767-006 
44767-007 
44767-008 
44767-009 
44767-010 
44767-011 
44767-012 
44767-013 
44767-014 
44767-015 
44767-016 
44767-017 
44767-018 
44767-019 
44767-020 
44767-021 
44767-022 
44767-023 
44767-024 
44767-025 
44767-026 
44767-027 
44767-028 
44767-029 
44767-030 
44767-031 
44767-032 
44767-033 
44767-034 
44767-035 
44767-036 
44767-037 
44767-038 
44767-039 
44767-040 
44767-041 
44767-042 
44767-043 
44767-044 
44767-045 
44767-046 
44767-047 
44767-048 
44768-001 
44768-002 
44768-003 
44768-004 
44768-005 
44768-006 
44768-007 
44768-008 
44768-009 
44768-010 
44768-011 
44768-012 
44768-013 


Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #1 
Wisconsin, USA (site #2 
Wisconsin, USA (site #2 
Wisconsin, USA (site #2 
Wisconsin, USA (site #2 
Wisconsin, USA (site #2 
Wisconsin, USA (site #2 
Wisconsin, USA (site #2 
Wisconsin, USA (site #2 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
Kentucky, USA 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 


Myotis septentrionalis 
Perimyotis subflavus 
Myotis septentrionalis 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Myotis septentrionalis 
Myotis septentrionalis 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Myotis septentrionalis 
Myotis lucifugus 
Perimyotis subflavus 
Myotis septentrionalis 
Myotis septentrionalis 
Myotis septentrionalis 
Myotis sp.* 

Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis sp.* 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Myotis grisescens 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Myotis sodalis 
Myotis sodalis 
Myotis lucifugus 
Perimyotis subflavus 
Myotis grisescens 
Perimyotis subflavus 
Myotis sodalis 
Myotis sodalis 
Myotis sodalis 
Myotis sodalis 
Myotis sodalis 
Myotis sodalis 
Myotis sodalis 
Myotis sodalis 
Myotis sodalis 
Myotis sodalis 
Myotis sodalis 
Myotis sodalis 
Perimyotis subflavus 
Perimyotis subflavus 
Perimyotis subflavus 
Myotis sodalis 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 


03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
03 March 2014 
10 March 2014 
10 March 2014 
10 March 2014 
10 March 2014 
10 March 2014 
10 March 2014 
10 March 2014 
10 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
04 March 2014 
19 March 2014 
19 March 2014 
19 March 2014 
19 March 2014 
19 March 2014 
19 March 2014 
19 March 2014 
19 March 2014 
19 March 2014 
19 March 2014 
19 March 2014 
19 March 2014 
19 March 2014 


no 
no 
no 
no 
no 
no 
yes 
yes 
no 
no 
no 
no 
no 
yes 
no 
no 
no 
no 
no 
no 
yes 
no 
no 
no 
no 
yes 
no 
yes 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
yes 
no 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
no 
yes 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
yes 
no 
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Table 1 (cont.) 


Individual identifier 


Location 


Host species 


Sampling date 


Malassezia vespertilionis isolated 


44768-014 
44768-015 
44768-016 
44768-017 
44768-018 
44768-019 
44768-020 
44768-021 
44768-022 
44768-023 
44768-024 
44768-025 
44768-026 
44768-027 
44768-028 
44768-029 
44768-030 
44768-031 
44769-001 
44769-002 
44769-003 
44769-004 
44769-005 
44769-006 
44769-007 
44769-008 
44769-009 
44769-010 
44769-011 
44769-012 
44769-013 
44769-014 
44769-015 
44769-016 
44769-017 
44769-018 
44769-019 
44769-020 
44769-021 
44769-022 
44769-023 
44769-024 
44769-025 
44797-032 
44797-033 
44797-034 
44797-035 
44797-036 
44797-037 
44797-038 
44797-039 
44797-040 
44797-041 
44797-042 
44797-043 
44797-044 
44797-045 
44797-046 
44797-047 
44797-048 
44797-049 
44797-050 
44797-051 
44797-052 
44797-053 
44797-054 
44797-055 
44797-056 
44797-057 
44797-058 
44797-059 
44797-060 
44797-061 
44797-062 
44797-063 
44797-064 
44797-065 
44797-066 
44797-067 
44797-068 
44797-069 
44797-070 
44797-071 
44797-072 
44797-073 
44797-074 
44797-075 
44797-076 
44797-077 


New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
New York, USA (site #1) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
Wisconsin, USA (site #3) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
New York, USA (site #2) 
Missouri, USA 

Missouri, USA 

Missouri, USA 

Missouri, USA 

Missouri, USA 

Missouri, USA 

Missouri, USA 

Missouri, USA 

Missouri, USA 

Missouri, USA 

Missouri, USA 

Missouri, USA 

Missouri, USA 

Missouri, USA 

Missouri, USA 

Missouri, USA 


Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
Myotis lucifugus 
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Table 1 (cont.) 


Individual identifier Location Host species Sampling date Malassezia vespertilionis isolated 
44797-078 Missouri, USA Perimyotis subflavus 24 February 2015 no 
44797-100 Wisconin, USA (site #4) Myotis septentrionalis 27 January 2015 yes 
44797-101 Wisconin, USA (site #4) Myotis septentrionalis 27 January 2015 no 
44797-102 Wisconin, USA (site #4) Myotis septentrionalis 27 January 2015 no 
44797-103 Wisconsin, USA (site #1) Myotis septentrionalis 28 January 2015 yes 
44797-104 Wisconsin, USA (site #1) Myotis septentrionalis 28 January 2015 yes 
44797-105 Wisconsin, USA (site #1) Myotis septentrionalis 28 January 2015 no 
44797-106 Wisconsin, USA (site #1) Myotis septentrionalis 28 January 2015 no 
44797-107 Wisconsin, USA (site #1) Myotis septentrionalis 28 January 2015 no 
44797-108 Wisconsin, USA (site #1) Myotis septentrionalis 28 January 2015 no 
44797-109 Wisconsin, USA (site #1) Myotis septentrionalis 28 January 2015 no 
44797-110 Wisconsin, USA (site #1) Myotis septentrionalis 28 January 2015 no 
44797-111 Wisconsin, USA (site #5) Myotis septentrionalis 29 January 2015 no 
44797-112 Wisconsin, USA (site #5) Myotis septentrionalis 29 January 2015 no 
44797-113 Wisconsin, USA (site #5) Myotis septentrionalis 29 January 2015 no 
44797-114 Wisconsin, USA (site #5) Myotis septentrionalis 29 January 2015 no 
44797-115 Wisconsin, USA (site #5) Myotis septentrionalis 29 January 2015 yes 
44797-116 Wisconsin, USA (site #5) Myotis septentrionalis 29 January 2015 no 
44797-123 Wisconsin, USA (site #6) Myotis septentrionalis 02 March 2015 no 
44797-124 Wisconsin, USA (site #6) Myotis septentrionalis 02 March 2015 yes 
44797-125 Wisconsin, USA (site #6) Myotis septentrionalis 02 March 2015 no 
44797-126 Wisconsin, USA (site #6) Myotis septentrionalis 02 March 2015 no 
44797-127 Wisconsin, USA (site #6) Myotis septentrionalis 02 March 2015 yes 
44797-128 Wisconsin, USA (site #6) Myotis septentrionalis 02 March 2015 no 
44797-129 Wisconsin, USA (site #6) Perimyotis subflavus 02 March 2015 no 
44797-130 Wisconsin, USA (site #6) Myotis septentrionalis 02 March 2015 no 
44797-131 Wisconsin, USA (site #6) Myotis septentrionalis 02 March 2015 yes 
44797-132 Alabama, USA Perimyotis subflavus 11 February 2015 no 
44797-133 Alabama, USA Myotis grisescens 11 February 2015 no 
44797-134 Alabama, USA Perimyotis subflavus 11 February 2015 no 
44797-135 Alabama, USA Myotis sodalis 11 February 2015 yes 
44797-136 Alabama, USA Myotis grisescens 11 February 2015 no 
44797-137 Alabama, USA Myotis grisescens 11 February 2015 yes 
44797-138 Alabama, USA Myotis grisescens 11 February 2015 yes 
44797-139 Alabama, USA Myotis grisescens 11 February 2015 yes 
44797-140 Alabama, USA Perimyotis subflavus 11 February 2015 no 
44797-141 Alabama, USA Myotis grisescens 11 February 2015 yes 
44797-142 Alabama, USA Myotis grisescens 11 February 2015 no 
44797-143 Alabama, USA Myotis grisescens 11 February 2015 yes 
44797-144 Alabama, USA Myotis grisescens 11 February 2015 no 
44797-145 Alabama, USA Perimyotis subflavus 11 February 2015 no 
44797-146 Alabama, USA Myotis grisescens 11 February 2015 no 
44797-147 Alabama, USA Myotis grisescens 11 February 2015 no 
44797-148 Alabama, USA Myotis grisescens 11 February 2015 yes 
44797-149 Alabama, USA Perimyotis subflavus 11 February 2015 no 
44797-150 Alabama, USA Myotis grisescens 11 February 2015 yes 
44797-151 Alabama, USA Myotis grisescens 11 February 2015 no 
44797-152 Alabama, USA Myotis grisescens 11 February 2015 yes 
44797-153 Alabama, USA Myotis sodalis 11 February 2015 yes 
44797-154 Alabama, USA Myotis grisescens 11 February 2015 no 
44797-155 Alabama, USA Myotis grisescens 11 February 2015 no 
45701-660 California, USA Myotis yumanensis 07 May 2016 no 
45701-661 California, USA Myotis sp.** 02 May 2016 no 
45701-663 California, USA Myotis yumanensis 18 April 2016 no 
45701-664 California, USA Myotis yumanensis 02 May 2016 yes 
45701-665 California, USA Myotis californicus 07 May 2016 yes 
45701-666 California, USA Myotis californicus 02 May 2016 no 
45701-668 California, USA Myotis yumanensis 09 May 2016 yes 
45701-669 California, USA Myotis yumanensis 18 April 2016 no 
45701-671 California, USA Myotis yumanensis 09 May 2016 no 
45701-672 California, USA Myotis californicus 02 May 2016 no 
45701-674 California, USA Myotis californicus 26 April 2016 no 
45701-675 California, USA Myotis yumanensis 07 May 2016 no 
45701-676 California, USA Myotis sp.** 02 May 2016 yes 
45701-677 California, USA Myotis yumanensis 18 April 2016 yes 
45701-678 California, USA Myotis californicus 18 April 2016 yes 
45701-680 California, USA Myotis sp.** 18 April 2016 no 
45701-681 California, USA Myotis sp.** 09 May 2016 no 
45701-682 California, USA Myotis californicus 19 April 2016 yes 
45701-683 California, USA Myotis yumanensis 20 April 2016 no 
45701-684 California, USA Myotis californicus 09 May 2016 no 
45701-685 California, USA Myotis yumanensis 20 April 2016 no 
45701-686 California, USA Myotis thysanodes 26 April 2016 yes 
45701-687 California, USA Myotis yumanensis 20 April 2016 no 
45701-688 California, USA Myotis yumanensis 20 April 2016 yes 
45701-689 California, USA Myotis sp.** 02 May 2016 no 
45701-691 California, USA Myotis californicus 16 May 2016 yes 
45701-696 California, USA Myotis californicus 16 May 2016 no 
45701-699 California, USA Lasionycteris noctivagans 16 May 2016 no 
45701-714 California, USA Lasionycteris noctivagans 16 May 2016 yes 
45701-719 California, USA Myotis californicus 16 May 2016 yes 
45704-161 Pennsylvania, USA Myotis leibii 26 February 2016 yes 
46375-001 California, USA Myotis californicus 08 May 2017 yes 
46375-002 California, USA Myotis californicus 08 May 2017 yes 
46375-003 California, USA Lasionycteris noctivagans 08 May 2017 yes 
46375-004 California, USA Lasionycteris noctivagans 08 May 2017 no 


* 


either Myotis lucifugus or Myotis septentrionalis. 


** either Myotis lucifugus or Myotis yumanensis. 
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Table 2 Fungal genomes and summary statistics used for whole genome phylogenic analysis. 


Species Strain GenBank Locus tag Assembly Number of Scaffold N50 (Kb) Percent GC Protein coding 
identifier accession number size (Mb) scaffolds gene models 
Malassezia caprae CBS 10434* GCA 001264625.1 MCA1 7.58 229 110 59.7896 3.553 
M. cuniculi CBS 11721* GCA 001264635.1 MCU1 7.459 76 522 58.99% 3.167 
M. dermatis CBS 9169* GCA_001264665.1 MDM1 7.54 111 189 59.10% 3.538 
JCM 11348 GCA_001600775.1 MDM2 7.551 18 1.325 58.98% 3.544 
M. equina CBS 9969* GCA_001264685.1 MEQ1 7.658 117 372 58.00% 3.232 
M. furfur JPLK 23 GCA_001265065.1 MFU6 7.79 2092 14 64.18% 3.087 
CBS 1878* GCA_001265055.1 MFU1 13.865 3460 14 63.91% 5.418 
CBS 4172 GCA_001264895.1 MFU2 14.347 3453 15 64.23% 5.672 
CBS 7019* GCA 001264875.1 MFU3 13.707 3262 15 63.7096 5.565 
CBS 7710 GCA 001264865.1 MFU4 15.232 4053 14 63.89% 5.835 
CBS 7982 GCA 001265045.1 MFU5 7.877 1694 20 64.06% 3.158 
M. globosa CBS 7966* GCA_001264805.1 MGL2 8.94 113 724 52.02% 4.245 
CBS 7874 GCA_001264815.1 MGL1 8.938 138 398 51.87% 4.191 
CBS 7990 GCA_001264795.1 MGL3 8.884 108 414 52.07% 3.703 
M. japonica CBS 9431* GCA_001264785.1 MJA1 8.341 295 66 62.38% 4.215 
JCM 11963 GCA_001600795.1 MJA2 8.364 16 814 62.33% 4.122 
M. nana CBS 9557* GCA_001265015.1 MNA1 7.607 95 492 57.93% 3.785 
JCM 12085 GCA_001600835.1 MNA2 7.579 13 1.323 57.96% 3.734 
M. obtusa CBS 7876* GCA_001264985.1 MOB1 7.842 1709 22 62.15% 2.893 
M. pachydermatis CBS 1879* GCA_001264975.1 MPA1 8.158 61 957 55.08% 4.134 
M. restricta CBS 7877* GCA_001264765.1 MRE1 7.249 90 402 55.83% 3.556 
CBS 8742 GCA_001264725.1 MRE2 7.26 69 666 55.79% 3.569 
M. slooffiae CBS 7956* GCA_001264965.1 MSL1 8.425 1641 15 65.82% 3.422 
M. sympodialis ATCC 42132 GCA 001264925.1 MSY1 7.546 824 54 58.77% 3.055 
ATCC 44340 GCA 001264715.1 MSY2 7.562 769 59 58.88% 3.080 
ATCC 96806 GCA_001264705.1 MSY3 7.526 1030 44 58.80% 3.946 
M. vespertilionis CBS 15041* GCA_002818225.1 MVES 7.581 14 844 56.62% 3.791 
M. yamatoensis CBS 9725* GCA 001264885.1 MYA1 8.106 49 1.447 49.6296 3.971 


* type or neotype strain. 


and trimmed using trimAl v. 1.4.rev15 (Capella-Gutiérrez et al. 
2009). A maximum-likelihood phylogeny was estimated using 
RAxML v. 8.2.10 (Stamatakis 2014) (PROTGAMMALG, 100 
bootstrap replicates). As a secondary method, a Bayesian 
phylogeny was estimated using MrBayes v. 3.2.6 (Ronquist 
& Huelsenbeck 2003) through the CIPRES Science Gateway 
(Miller et al. 2010). For the Bayesian analysis, an LG model 
with gamma-distributed rate variation across sites was used. To 
generate the 50 % majority rule consensus tree, two runs, each 
with 1 000 000 generations and four chains, were performed. 
The chains were sampled every 250 generations with the first 
25 96 of sampled values discarded as burn-in. 


Multilocus sequence analysis 


To determine whether the isolates of Malassezia from bats rep- 
resented a single species, seven loci (in addition to ITS) from 
12 isolates (Table 3) were also examined: the D1/D2 region of 
26S rDNA (hereafter referred to as D1/D2), the 18S rDNA, chitin 
synthase CHS2, second largest subunit of RNA polymerase II 
(RPB2), B-tubulin (B-tub), translation elongation factor EF-1a 
(TEF1), and minichromosome maintenance complex compo- 
nent 7 (MCM7). DNA was extracted using the methods of Lorch 
etal. (2015). The D1/D2 region was amplified using primers NL4 
(5-GGT CCG TGT TTC AAG ACG G-3) and NL1 (5-GCA TAT 
CAA TAA GCG GAG GAA AAG-3’) (O'Donnell 1993); cycling 
conditions: 94 °C for 5 min; 30 cycles of 94 °C for 45 s, 51 °C 
for 1 min, and 72 °C for 3 min; and a final extension of 72 °C 
for 10 min. The 18S rDNA was amplified with forward (5’-ATC 
TGGTTGATC CTG CCAGT-3’) and reverse (5-TCC TCC GCT 
TAT TGA TAT GC-3) primers described by Sugita & Nakase 
(1999); cycling conditions: 94 °C for 3 min; 30 cycles of 94 °C 
for 30 s, 55 °C for 1 min, and 72 °C for 2 min and 30 s; and 
a final extension of 72 °C for 10 min. A portion of CHS2 was 
amplified with primers ChiSyn2f (5-CTG AAG CTT CAN ATG 
TAY AAY GAR GAY-3’) and ChiSyn2r (5--GTT CTC GAG YTT 
RTAYTC RAA RTT YTG-3’) (Bowen et al. 1992, Cabaries et al. 


2007); cycling conditions: 94 °C for 5 min; 45 cycles of 94 °C for 
1 min, 50 °C for 2 min, and 72 °C for 3 min; and a final exten- 
sion of 72 °C for 6 min. A fragment of RPB2 was amplified with 
primers fRPB2-5F (5'-GAY GAY MGW GAT CAY TTY GG-3’) 
and fRPB2-7cR (5'-CCC ATR GCT TGY TTR CCC AT-3’) (Liu 
et al. 1999); cycling conditions: 94 °C for 4 min; 40 cycles of 
94 °C for 1 min, 50 °C for 1 min, and 72 °C for 1 min; and a final 
extension of 72 °C for 8 min. A portion of B-tub was amplified 
with primers F-Btub (5-CAR GCY GGT CAR TGY GGT AAC 
CA-3’) and F- Btub4r (5-GCC TCA GTR AAY TCC ATY TCR 
TCC AT-3’) (Einax & Voigt 2003); cycling conditions: 95 °C for 
5 min; 30 cycles of 95 °C for 30 s, 50 °C for 1 min, and 72 °C 
for 1 min; and a final extension of 72 °C for 10 min. A portion of 
TEF1 was amplified with primers EF1-983F (5'-GCY CCY GGH 
CAY CGT GAY TTY AT-3’) and EF1-2218R (5-ATG ACACCR 
ACR GCR ACR GTY TG-3’) (Rehner & Buckley 2005); cycling 
conditions: 94 °C for 2 min; 47 cycles of 94 °C for 1 min, 53 °C 
for 1 min, and 72 °C for 1 min and 40 s; and a final extension of 
72 °C for 10 min. A portion of MCM7 was amplified with primers 
MCM7-709 (5'-ACN MGN GTN TCV GAY GTH AAR CC-3’) 
and MCM7-1348 (5”-GAY TTD GCN ACN CCN GGR TCW 
CCC AT-3’) (modified slightly from Schmitt et al. 2009); cycling 
conditions as described for TEF1. Reactions were carried out 
using GoTaq® Flexi DNA polymerase (Promega Corporation, 
Madison, WI) according to the manufacturer's instructions (with 
final concentrations of 1.5 mM MgCl, solution and 1 uM of each 
primer) except that twice the recommended amount of Taq poly- 
merase (2.5 U) and 1—3 uL of template were added per 50 uL 
reaction. When necessary, PCR products were gel-purified prior 
to sequencing using the QlAquick Gel Extraction Kit (Qiagen 
Inc., Valencia, CA). All PCR products were sequenced in both 
directions using the same primers described for the initial 
amplifications. Additional internal sequencing primers were 
used for some loci: 18S rDNA, forward primers (5-GCT ACC 
ACA TCC AAG GAA GG-3’, 5-CTG CGA AAG CAT TTG CCA 
AGG-3’, 5-TCT GGG CCG CAC GCG CGC TAC ACT G-3’) 
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Table3 (cont.) 


GenBank accession numbers 


p-tub 


RPB2 MCM7 CHS2 


TEF1 


18S rDNA 


ITS D1/D2 


Species 


Other identifier 


Isolate identifier 


NA NA NA NA NA NA 


NA 


MF669457 


M. vespertilion. 
M. vespertilion. 
M. vespertilion. 
M. vespertilion: 


NVVHC 45701-664 


MF669382 MF669334 MF669322 MF669358 MF669370 MF669346 


MF669394 


NA 
NA 
NA 
NA 


MF669458 


is 
is 


CBS 15048 


NWHC 45701-665 


MF669459 
MF669460 
MF669461 
MF669462 


is 
is 


M. vespertilionis 


M. vespertilionis 


NWHC 45701-668 
NWHC 45701-676 
NWHC 45701-677 
NWHC 45701-678 


MF669383 MF669335 MF669323 MF669359 MF669371 MF669347 


MF669395 


MF669463 


CBS 15051 
CBS 15049 


NWHC 45701-682 


MF669348 
NA 


MF669372 
NA 


MF669360 
NA 


MF669324 
NA 


MF669336 
NA 


MF669384 
NA 


MF669396 


MF669464 


is 
is 
is 
is 


NWHC 45701-686 


MF669465 


M. vesperti 
M. vespertil 
M. vesperti 
M. vespertil 


NVVHC 45701-688 


MF669385 MF669337 MF669325 MF669361 MF669373 MF669349 


MF669397 


NA 
NA 
NA 
NA 
NA 


MF669466 


on: 
on: 


CBS 15050 


NVVHC 45701-691 


MF669467 
MF669468 
MF669469 
MF669470 
MF669471 


S 
S 


NWHC 45701-714 
NWHC 45701-719 
NWHC 45704-161 
NWHC 46375-001 
NWHC 46375-002 
NWHC 46375-003 
M9927* 


MF669386 MF669338 MF669326 MF669362 MF669374 MF669350 


MF669398 
AB070361 


MF669472 
AB070356 


CBS 15047 


M. vespertilion 
M. dermatis 
M. furfur 

M. furfur 

M. furfur 


LFFX01000059.1 
LFGIO1003290.1 


LFFX01000013.1 LFFX01000105.1 LFFX01000023.1 
LFGIO1002103.1 


LFGIO1000568.1 


LFFX01000095.1 
LFGIO1000839. 1 


KF706452.1 


CBS 9169* 


LFGG01000889.1 


N 


LFGG01003117.1 


LFGIO1001900.1 
LFGG01002621.1 


LFGG01000118.1 


LFGG01001873.1 


EU192363.1 


AY743603 NA 


AY743602 


AY743634 
AY743635 


CBS 1878* 
CBS 7019* 


A 


AY083223.1 


NCPF 3349 
CBS 9431* 


LFGC01000917.1 
LFGB01000046.1 


LFDB01000111.1 
NA 


LFGC01001067.1 
LFGB01000036.1 


LFDB01000119.1 
NA 


LFDB01000183.1 
LFGC01000795.1 


NA 
LFGB01000009.1 


LFDB01000004.1 
LFGC01000155.1 


NA 
LFGB01000018.1 


LFDBO01000216.1 
LFGC01001328.1 
LFGB01000029.1 


KF706458.1 
EU192365.1 
LFGB01000057.1 


NA 


EF140672 
AY743629 


NA 
AY743605 


EF 140669 
AY387137 


NA 
AY387139 


M. pachydermatis 


M. japonica 
M. obtusa 
M. obtusa 


CBS 7876* 
CBS 7968 
CBS 1879* 


sequence data not available or not used for analyses. 


NA= 


* type or neotype strain. 


and reverse primers (5-TGG AAT TAC CGC GGC TGC TGG 
CAC C-3’, 5-TCC TTG GCAAAT GCT TTC GCA G-3’, 5-CCG 
TCAATT CCT TTAAGT TTC AGC C-3’, 5-AAG GTC TCG TTC 
GTT ATC G-3’, 5'-GAC GGG CGG TGT GTA CAA AGG GCA 
G-3’) (Sugita 8, Nakase 1999), RPB2, RPB2-6f (5-TGG GGH 
ATG GTV TGY CCB GC-3’) and RPB2-6r (5-GCV GGR CAB 
ACC ATD CCC CA-3”) (modified slightly from Liu et al. 1999); 
and TEF1, EF1-1577F (5'-CAG GAY GTN TAC AAG ATY GGT 
GG-3”) and EF1-1567R (5-ACH GTR CCRATA CCA CCR ATC 
TT-3’) (Rehner & Buckley 2005). 


A phylogenetic analysis was conducted using newly generated 
sequences from the bat-associated isolates of Malassezia and 
from existing sequence data in GenBank for type cultures of 
a subset of the recognised species of Malassezia for which 
sufficient sequence data were available (Table 3). Members 
of the genus residing in the same core clade as the type iso- 
late (as determined by the whole-genome analysis described 
above) were included in this analysis as were representatives 
from the other two core clades described by Wu et al. (2015). 
Sequence data for protein-coding genes were obtained from 
whole genome sequences deposited in GenBank by Wu et 
al. (2015), while that for multicopy genes originated from vari- 
ous sources (see Table 3). For M. obtusa and one isolate of 
M. furfur, 18S sequence data were not available for the type 
isolates; instead, sequence data from non ex-type strains were 
substituted for that locus. 


Nucleotide sequences were aligned independently for each 
locus using MUSCLE in MEGA v. 6 (Tamura et al. 2013), and 
all gaps were deleted. MEGA 6 was also used to determine the 
best substitution model for each locus. A multigene phylogenetic 
analysis was then conducted by concatenating the final align- 
ments of all eight loci. A Bayesian analysis was run as described 
above, except that 5000000 generations were used for each 
run and the sampling frequency was set to 1000. Data was 
partitioned by locus and (for coding genes) by nucleotide codon 
position. A Kimura 2-parameter model with gamma distribution 
was applied to the non-coding loci (ITS, 18S rDNA, and D1/D2); 
a Kimura 2-parameter model with gamma distribution and invari- 
ant sites was used for CHS2; a general time-reversible model 
with gamma distribution was used for B-tub, TEF1, and RPB2; 
and an HKY model with gamma distribution and invariant sites 
was applied to MCMT. 


Physiological and morphological 
characterisation of isolates 


Seven isolates of the bat-associated Malassezia sp. that were 
analysed genetically were also characterised physiologically 
and morphologically using criteria commonly employed to dis- 
tinguish species in the genus (Guého-Kellerman et al. 2010). To 
determine the influence of temperature on growth, all isolates 
were incubated at various temperatures as described by Guého 
et al. (1996). Due to the fastidious nature ofthe isolates, growth 
temperature experiments were conducted on LNA instead of 
modified Dixon's agar (mDA; Guillot et al. 1996). Growth was 
assessed at the following temperatures: 7, 24, 30, 37, and 
40 °C. Inclusion of growth characteristics at 7 °C and 24 °C 
are not standard for Malassezia, but were performed because 
the fungus was isolated from hibernating bats, which maintain 
low body temperatures. Plates were inoculated by transferring 
cells with an inoculating loop and streaking for isolated colonies. 
The diameter of isolated colonies was measured and colony 
morphology descriptions were recorded every 10 d for a total 
of 50 d. Cell morphology was assessed from wet mounts with 
lactophenol cotton blue stain conducted on 10-d-old cultures 
that were grown on LNA at 24 °C. 


The ability of the bat-associated Malassezia isolates to grow on 
mDA and without lipid supplementation on Sabouraud dextrose 
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agar (SDA) was tested by inoculating these media as described 
above. Utilisation of different types of Tween (i.e., 0.1 % Tween 
80, 0.5 % Tween 40, 0.5 % Tween 60, and 10 % Tween 20) was 
tested according to the methodologies of Guého et al. (1996). 
The seven isolates were also characterised using the Tween 
Diffusion and Cremophor EL Assimilation Tests (Guillot et al. 
1996). For all of these experiments, cultures were incubated at 
32 °C and checked every seven days for a total of 50 d. Tests 
were also performed at 24 *C to ensure that lack of positive 
results was not due to incubation at temperatures outside the 
optimal growth range for the Malassezia sp. 


Additional physiological tests included catalase reaction and 
P-glucosidase activity. The catalase reaction was performed by 
harvesting a loop full of cells from a 7-d-old culture grown on 
LNAat 24 °C, smearing the cells onto a glass slide, and adding 
one drop of 3 % hydrogen peroxide (Guého et al. 1996). The 
ability to hydrolyse esculin (i.e., B-glucosidase activity) was 
tested following the methods of Mayser et al. (1997) except 
that 7-d-old cultures grown on LNA were used to inoculate the 
medium due to the slower growth rate of the bat-associated iso- 
lates relative to other species of Malassezia. The esculin agar 
tubes were incubated at 32 °C and checked daily for 14 d, and 
thereafter they were checked every seven days for an additional 
30 d. The test was considered positive if a black precipitate was 
produced. Mating experiments were not performed. 


RESULTS 


Isolation of Malassezia from bats 


Fungal colonies that resembled Malassezia were visible on LNA 
medium after 40—50 d of incubation at 7 °C. These putative 
Malassezia colonies were observed growing in culture from sam- 
ples collected from 75 of the 264 (28 96) bats examined. These 
included nine host species from seven U.S. states (Table 1; 
Fig. 2). The ITS sequences of 74 isolates shared 99.7—100 % 
sequence identity with one another across the approximately 


760-bp fragment that was analysed. The remaining isolate had 
an ITS sequence divergent from the other 74 isolates and was 
not included in further analyses. When searched against the 
GenBank database, ITS sequences of the 74 bat-associated 
isolates most closely matched M. japonica (92 % sequence 
identity) within the portion of the sequences corresponding 
to the 5.8S rDNA; however, ITS1 and ITS2 were highly di- 
vergent from sequences that resided within GenBank. The 
ITS sequence data are available in GenBank (MF669451— 
MF669472). 


Some isolates of the Malassezia sp. lost vigour and eventually 
failed to grow after several passes on LNA. Transfer of these 
isolates to mDA did not cause reinvigoration but rather seemed 
to facilitate further decline of the cultures. The reason for this 
was not determined. Despite losing some isolates, the major- 
ity grew well in the laboratory. Representative isolates were 
deposited in the Westerdijk Fungal Biodiversity Institute and 
the UAMH Centre for Global Microfungal Biodiversity culture 
collections (Table 3). 


Whole genome sequence analysis 


The Whole Genome Shotgun project was deposited at DDBJ/ 
ENA/GenBank under the WGS Project PECA00000000 (ac- 
cession GCA 002818225.1). The version described in this 
paper is version PECA01000000. Raw sequencing data is 
available from the NCBI SRA via the accession SRP121079. 
The genome of the Malassezia isolated from bats was 7.581 
Mb, contained in 14 scaffolds, which is consistent with other 
Malassezia species (Boekhout et al. 1998, Wu et al. 2015). The 
annotated genome is estimated to contain 3791 protein coding 
gene models (Table 2). Twenty-seven genomes of Malassezia 
species were obtained from NCBI GenBank, however most of 
these data contained only assemblies without annotated gene 
models. Thus, all genome assemblies were re-annotated us- 
ing funannotate in order to generate comparable annotations 
between species. The number of predicted genes was similar 
to previously published reports (Table 2) (Wu et al. 2015). 
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Fig. 2 Map of the United States, showing the locations of sampling sites from which bats yielded isolates of Malassezia vespertilionis sp. nov. States from 
which isolates were obtained are labelled (AL = Alabama, CA = California; KY = Kentucky; MO = Missouri; NY = New York; PA = Pennsylvania; WI = Wisconsin). 
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1/100 M. slooffiae CBS 7956 
M. cuniculi CBS 11721 
M. pachydermatis CBS 1879 
1/96 1/100, M. nana CBS 9557 
M. nana JCM 12085 
ə M. equina CBS 9969 
M. caprae CBS 10434 
1/100 M. dermatis CBS 9169 
4/100 ə M. dermatis JCM 11348 
1/100|., ; M. sympodialis ATCC 42132 Clade B 
100 M. sympodialis ATCC 44340 
1/100] M sympodialis ATCC 96806 
1/100, M. restricta CBS 7877 
4/100 M. restricta CBS 8742 
1/100 17100 M. globosa CBS 7874 
M. globosa CBS 7966 
1/100 L M. globosa CBS 7990 
M. vespertilionis sp. nov. CBS 15041 
1/100 | M. japonica CBS 9431 
1/100 M. japonica JCM 11963 
iog M. obtusa CBS 7876 
M. yamatoensis MY9725 
11967 M. furfur CBS 1878 Clade A 
196- M. furfur CBS 4172 


1/100 M. furfur CBS 7019 
7 "M. f rfur CBS 7710 
anoo] “runur 


M. furfur CBS 7982 
M. furfur JPLK23 


Clade C 


1/100 


1/- 


Ustilago maydis 521 


0.09 


Fig. 3 Phylogenetic tree of the genus Malassezia based on concatenated amino acid sequences of 254 conserved orthologues. The tree from the Bayesian 
analysis is shown, but the tree generated from the maximum likelihood analysis had an identical topology. Posterior probabilities (Bayesian)/bootstrap values 
(maximum likelihood), respectively, are shown at the nodes. Ustilago maydis was used to root the tree. Clades A, B, C as described by Wu et al. (2015) are 


illustrated. Based on the analyses, M. vespertilionis sp. nov. is a basal member of clade A. 


1 M. pachydermatis CBS 1879 


CBS 15042 
CBS 15043 


CBS 15046 
CBS 15044 
CBS 15045 
CBS 15041 
CBS 15047 
CBS 15050 
CBS 15048 
CBS 15051 
CBS 15049 
M. japonica CBS 9431 

M. obtusa CBS 7876* 


11M. furfur CBS 1878 
M. furfur CBS 7019 
M. cuniculi CBS 11721 


M. dermatis CBS 9169 


UAMH 11925 


M. vespertilionis sp. nov. 


M. yamatoensis CBS 9725 


0.03 


Fig.4 Phylogenetic tree resulting from a Bayesian analysis of concatenated nucleotide sequences from eight loci (ITS, 18S rDNA, D1/D2 region, and portions 
ofthe B-tub, TEF1, MCM7, RPB2, CHS2 genes) of 12 Malassezia isolates from bats, all Malassezia species from clade A for which sufficient genetic data was 
available, and representative members from clades B and C (Wu et al. 2015). Posterior probabilities are presented at each node. All examined isolates from 
bats formed a well-supported clade, suggesting that they represent a single taxon referred to herein as M. vespertilionis sp. nov. 
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Using the annotated genomes and Ustilago maydis as an 
outgroup, we identified core conserved gene models in each 
genome using BUSCO v. 2 orthologous groups and generated 
a concatenated alignment of 254 gene models found in all 29 
genomes studied. Using both maximum likelihood and Bayes- 
ian methods, the bat-associated Malassezia was placed as the 
most basal member of clade A (Fig. 3). 


Multilocus sequence analysis 


To demonstrate that all of the Malassezia isolates from bats 
represented a single taxon, multilocus sequencing analysis 
was performed on 12 isolates (including the type isolate). The 
isolates for the analysis were selected to represent a broad 
range of host species, geographic locations, and strains with 
slight sequence variations in the ITS region. The portions se- 
quenced of D1/D2 region, B-tub, TEF1, and RPB2 were 100 % 
identical among the 12 isolates, and the ITS region, 18S rDNA, 
MCMT, and CHS2 shared at least 99.7, 99.9, 99.8, and 99.8 % 
sequence identity, respectively, between the isolates examined. 
The sequences generated for the multilocus analysis are avail- 
able in GenBank (D1/D2: MF669394 — MF669398; 18S rDNA: 
MF669382-MF669386, B-tub: MF669334 - MF669338; TEF1: 
MF669322-MF669326, RPB2: MF669358 - MF669362; MCMT: 
MF669370—MF669374; CHS2: MF669346— MF669350). 


The final alignments used for the phylogenetic analysis con- 
sisted of the following numbers of characters: ITS region, 402 
characters; 18S rDNA, 1065 characters; D1/D2 region, 541 
characters; B-tub, 1041 characters; TEF1, 987 characters; 
MCM7, 600 characters; RPB2, 1086 characters; and CHS2, 
534 characters. The small number of characters included in 
the ITS region alignment was due to high divergence between 
the different Malassezia species within the ITS1 and ITS2 
regions, resulting in their subsequent removal from the align- 
ment. Thus, the final 'ITS' alignment consisted almost entirely 
of sequence data representing the 5.8S rDNA which is more 
highly conserved among species of Malassezia. The Bayesian 
analysis based on the eight concatenated sequences produced 
a tree that showed the same relationships among the subset 
of species included as the phylogenetic analysis based on 
whole genome sequencing. The 12 Malassezia isolates from 
bats all grouped together into a well-supported clade that likely 
represented a single species (Fig. 4). 


Physiological and morphological 

characterisation of isolates 

The Malassezia sp. isolated from bats grew at all temperatures 
tested (7, 24, 30, 37, and 40 °C), with best growth (i.e., largest 
colony diameters) occurring at 24 °C. At higher temperatures, 
growth was slower. This was problematic since standard phy- 


siological tests used to characterise species of Malassezia are 
typically conducted at 32 °C for incubation periods that were too 
brief to allow for sufficient growth of this novel taxon from bats 
(Guého et al. 1996). Thus, we conducted most tests at both 
24 °C and 32 °C and allowed cultures to incubate for up to 50 d. 
Test results at the two temperatures were generally equivalent, 
although positive test results often required longer incubation 
times and produced weaker results at 32 °C compared to 24 °C. 
Detailed descriptions of growth at different temperatures, cell 
morphology, and colony morphology are presented in the spe- 
cies description and in Table 4. Colony and cell morphology 
are also shown in Fig. 5. 


The bat-associated Malassezia grew on mDA, but some isolates 
began to lose vigour (even more so than previously described 
on LNA) after several transfers on the medium. All isolates 
were lipid-dependent, failing to grow on SDA. Growth occurred 
on a variety of tween lipid sources, but not in the presence of 
Cremophor EL. The isolates were catalase and B-glucosidase 
negative. More detailed information is provided in the species 
description and in Table 4. 


SPECIES DESCRIPTION 


Malassezia vespertilionis J.M. Lorch & Vanderwolf, sp. nov. — 
MycoBank MB822382; Fig. 5 


Etymology. The species epithet refers to the host from which the fungus 
was isolated (n. vespertilio, Latin for bat; gen. n. vespertilionis, of a bat). 


Holotype. USA, Wisconsin, swab of wing skin of hibernating Myotis 
septentrionalis, 28 Jan. 2014, J.P. White (U.S. National Fungus Collections 
BPI 910536; culture ex-type CBS 15041 = UAMH 11924). 


Colonies are approximately 0.5—1.0 mm diam after 10 d of 
growth at 24 °C on LNA; 2-5 mm diam after 40 d. At 10 d, 
colonies are cream-coloured, flatto slightly convex, somewhat 
glossy, have entire margins, and have a crumbly or waxy con- 
sistency. At 40 d, colonies have irregular margins, are slightly 
raised, and have a prominent papilla near the centre (Fig. 5). 
Cells are ellipsoid or ovoid to (rarely) globose, ranging in size 
from 2—3 x 2-4 um (typically 2 x 3 um) (Fig. 5). Buds are 
formed monopolarly, usually on a narrow base. Growth occurs 
(sometimes poorly) on mDA; no growth observed on SDA. 
Isolates are catalase and B-glucosidase negative. Variable in 
lipid utilisation: no growth observed for Cremophor EL, weak 
or no growth for Tween 20, usually weak growth for Tween 80, 
weak to good growth for Tween 40, and good growth for Tween 
60. Growth occurs across a range of temperatures on LNA, but 
is slower at temperatures above and below 24 °C; specifically, 
individual pinpoint colonies are visible by day 30 at 7 °C and 
30 °C, and on day 40 at 37 °C. Growth is evident at 40 d for 


Fig. 5 Colony and cell morphology of M. vespertilionis sp. nov. grown on Leeming and Notman Agar at 24 °C. a. Colony size and morphology after 10 d of 
growth; b. colony size and morphology after 40 d of growth; c. cell morphology of 10-d-old culture. — Scale bars: a, b = 4 mm; c = 5 um. 
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cultures grown at 40 °C on areas of the medium where cells are 
placed at high densities; however, individual colonies are not 
grossly discernible. A sexual state was not observed; however, 
mating experiments were not explicitly performed. 


Additional specimens examined for physiological characteristics and multi- 
locus DNA sequencing. USA, Wisconsin, swab of wing skin of a hibernating 
Myotis septentrionalis, 3 Mar. 2014, J.P. White, CBS 15042; Kentucky, swab 
of wing skin of hibernating Myotis sodalis, 4 Mar. 2014, M.L. Verant, CBS 
15043; Kentucky, swab of wing skin of hibernating Myotis sodalis, 4 Mar. 
2014, M.L. Verant, UAMH 11925; New York, swab of wing skin of hibernating 
Myotis lucifugus, 19 Mar. 2014, M.L. Verant, CBS 15044; New York, swab of 
wing skin of hibernating Myotis lucifugus, 15 Jan. 2015, M.L. Verant, CBS 
15045; Alabama, swab of wing skin of hibernating Myotis grisescens, 11 Feb. 
2015, N. Sharp, CBS 15046. 

Additional isolates for which multilocus DNA sequencing was conducted. 
USA, California, swab of wing skin of Lasionycteris noctivagans, 8 May 2017, 
T.J. Weller, CBS 15047; California, swab of wing skin of Myotis californicus, 
7 May 2016, T.J. Weller, CBS 15048; California, swab of wing skin of Myotis 
thysanodes, 26 Apr. 2016, T.J. Weller, CBS 15049; California, swab of wing 
skin of Myotis californicus, 16 May 2016, T.J. Weller, CBS 15050; California, 
swab of wing skin of Myotis californicus, 19 Apr. 2016, T.J. Weller, CBS 15051. 


DISCUSSION 


The diversity of species within the genus Malassezia has been 
expanded in recent years due to increased sampling and use 
of genetic and molecular tools for distinguishing taxa. In the 
current study, we report on a species of Malassezia that is 
relatively common on skin of North American bats (i.e., cul- 
tured from 28 % of 264 individuals sampled). Isolates of this 
Malassezia sp. from bats were physiologically and genetically 
similar to one another. Specifically, sequences of four loci (D1/ 
D2 region, B-tub, TEF1, and RPB2) were identical between 
isolates. Minimal variation in the ITS region (99.7 % sequence 
identity), 18S rDNA (99.9 % sequence identity), and portions 
of the MCM7 (99.8 % sequence identity) and CHS2 (99.8 % 
sequence identity) genes were within the intraspecific variation 
documented in other species of Malassezia (e.g., Makimura et 
al. 2000, Hirai et al. 2004, Cabafies et al. 2007). As a group, 
the bat-associated isolates were highly similar to one another, 
and they were sufficiently distinct from all other known taxa. 
Thus, we propose that these bat-associated isolates represent 
the novel species M. vespertilionis. 


All isolates of M. vespertilionis subjected to physiological tests 
were catalase negative, a trait shared only with M. arunalokei, 
M. restricta, and some strains of M. furfur and M. pachydermatis 
(Guého et al. 1996, Guillot et al. 1998, reviewed by Batra et 
al. 2005, Honnavar et al. 2016). In contrast to M. arunalokei 
and M. restricta, M. vespertilionis is capable of growth (albeit 
slow) at 40 *C and can utilise multiple types of Tween. It can 
be distinguished from catalase negative strains of M. pachy- 
dermatis based on ability of the latter to grow on SDA without 
the supplementation of lipids. Differentiation of M. vespertili- 
onis from M. furfur based on physiological tests alone may 
be problematic due to reported variation in M. furfur (Batra et 
al. 2005). We examined a number of isolates of M. vespertili- 
onis and found some variability in results of the physiological 
tests for this species as well. Interpretation of physiological 
tests can be challenging (Gupta et al. 2004), and the inability 
of M. vespertilionis to grow sufficiently to produce positive 
results under the standard incubation procedures set forth by 
Guého et al. (1996) further complicates the use of these tests 
for identification. Thus, although physiological tests can often 
separate M. vespertilionis from other described taxa and may 
be helpful for some applications, we encourage the use of DNA 
sequencing (e.g., ITS region) to confirm identification. 


To date, Malassezia species have been isolated primarily from 
euthermic animals that maintain constant core body tempera- 


ture near 37 °C. This is consistent with the statement made by 
Guého-Kellerman et al. (2010) that members of the genus do 
not endure temperatures below 28 *C. Most of the bat species 
from which M. vespertilionis was isolated hibernate for up to 
seven months of the year at which time their body temperature 
is close to that of the surrounding environment. Specifically, 
Myotis lucifugus, Myotis sodalis, and Myotis septentrionalis 
were reported to prefer winter hibernacula with average air 
temperatures around 7.2, 8.5, and 9.1 °C, respectively (Brack 
2007). In the active season, a bat's body temperature may 
exceed 40 *C; however, bats often use bouts of torpor even 
during the active season, such that their body temperatures fre- 
quently fluctuate, sometimes approaching ambient temperature 
(Hock 1951, Studier 1981, Willis & Cooper 2009). Thus, the 
skin temperature of bats is highly variable, and this may explain 
the wide range of temperatures under which M. vespertilionis 
is capable of growth (from at least 7 °C to 40 °C). The ability 
to grow at such cool temperatures is noteworthy and may be 
unique to M. vespertilionis among the Malassezia. However, the 
lower growth limits of most Malassezia species have not been 
expressly described in literature, making comparisons difficult. 
The detection of DNA of uncultured malassezia-like organisms 
on corals and in terrestrial and marine environments suggests 
that other undescribed species of Malassezia may be capable 
of growth at ambient temperatures or reside on poikilothermic 
hosts (reviewed by Amend 2014). 


Previous phylogenetic studies of Malassezia have demon- 
strated disparities in the relationships between species when 
different loci were analysed (Cabafies et al. 2007, Castellá et 
al. 2014). This genealogical discordance has made decipher- 
ing taxonomic relationship among all members of the genus 
difficult. Wu et al. (2015) was able to better resolve Malassezia 
phylogeny through whole genome sequencing and concatena- 
tion of 164 genes. This phylogeny resulted in three main clades: 
clade A consists of M. furfur, M. obtusa, M. yamatoensis, and 
M. japonica; clade B contains M. sympodialis, M. dermatis, 
M. caprae, M. equina, M. nana, M. pachydermatis, M. globosa, 
and M. restricta; and clade C is comprised of M. cuniculi and 
M. slooffiae. Using the sequence data generated by Wu et al. 
(2015) and the whole genome sequence of M. vespertilionis 
produced in this study, we conducted a phylogenetic analysis 
using amino acid sequences of 254 core genes. The result- 
ing phylogenetic tree was identical to that of Wu et al. (2015) 
except that our analysis suggested that M. obtusa is basal 
to the subclade consisting of M. furfur and M. yamatoensis 
(Fig. 3). The Malassezia tree of Wu et al. (2015) had placed 
M. yamatoensis basal to M. furfur and M. obtusa; however, in 
that study the node was less well-supported compared to other 
species-level relationships. The slight change in tree topology 
and greater support for relationships within clade A in our tree 
may be due to the inclusion of more genes and M. vespertilionis 
in the analysis. These findings suggest that while sequencing 
loci traditionally used to differentiate species of Malassezia 
(e.g., ITS, D1/D2, B-tub, and CHS2 (Cabafies et al. 2007, 
2011, 2016, Castellá et al. 2014)) can be useful in identifying 
novel species, whole genome sequencing may be necessary to 
generate enough genetic data to fully resolve the relationship 
of those novel species with existing taxa. 


Our phylogenetic analyses indicate that M. vespertilionis is 
the most basal member of clade A. Malassezia arunalokei, 
M. brasiliensis, and M. psittaci were not included in the analyses 
because genetic data is available for only three loci each for 
these newly-described species (Cabafies et al. 2016, Honnavar 
et al. 2016). Based on existing sequence data, M. arunalokei, 
M. brasiliensis, and M. psittaci share only about 75—87 %, 
90—91 %, and 78—84 % DNA sequence identity, respectively, 
with M. vespertilionis in the ITS region, D1/D2 region, and 
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portion of the B-tubulin gene. Furthermore, previous analyses 
indicate that M. arunalokei is a member of clade B (Honnavar 
et al. 2016). Malassezia brasiliensis and M. psittaci are sister 
taxa to M. furfur and M. yamatoensis, respectively, which are 
both divergent from M. vespertilionis (Cabafes et al. 2016). 


The genus Malassezia may be more diverse than currently 
documented due to the difficulty of transporting and culturing 
many fragile and fastidious members of the genus, the historic 
use of morphological and physiological characteristics as the 
sole criteria to identify species (which can fail to distinguish 
cryptic species), and a lack of sampling of diverse taxonomic 
host groups (Amend 2014, Cabafes 2014). In the only other 
published study in which bats were specifically surveyed for 
Malassezia, Gandra et al. (2008) reportedly cultured M. furfur, 
M. globosa, M. pachydermatis, and M. sympodialis (based on 
physiological and morphological characteristics) from Pallas' 
mastiff bats (Molossus molossus) in Brazil. No DNA sequence 
data were generated and the isolates were apparently not 
deposited in a public culture collection, making it difficult to 
ascertain their true species assignments. However, because 
all isolates from Gandra et al. (2008) were either catalase posi- 
tive or lipid-independent (representing significant physiological 
deviations from M. vespertilionis) it may be that M. vespertilionis 
has not previously been isolated due to a lack of sampling effort 
of temperate bat species. Indeed, few studies have examined 
the fungal communities associated with bats, and those that 
have did not utilise fungal growth media suitable forthe isolation 
of lipid-dependent species of Malassezia (Grose & Marinkelle 
1966, Larcher et al. 2003, Voyron et al. 2011, Johnson et al. 
2013, Vanderwolf et al. 2013). Njus (2014) detected Malasse- 
zia spp. on the skin of bats by conducting fungal community 
analyses with next generation sequencing. Based on reported 
sequence identities with other species of Malassezia, at least 
some of the Malassezia detected by Njus (2014) may repre- 
sent M. vespertilionis. However, DNA sequence data from that 
project were not available in GenBank at the time of our study 
to confirm this. 


Malassezia yeasts are best known for their association with 
certain skin ailments. For example, M. pachydermatis has been 
implicated as the cause of dermatitis in (among other species) 
rhinoceroses (Bauwens et al. 1996), dogs (Gustafson 1955, 
Bond et al. 2004), and sea lions (Guillot et al. 1998, Nakagaki 
et al. 2000). More often, the role of Malassezia in skin diseases 
of animals is unknown, and targeting skin lesions for culture 
and PCR analyses likely masks the frequency with which these 
fungi act as commensals. For example, Neves et al. (2017) 
isolated Malassezia from a relatively high percentage (32.8 96) 
of free-ranging golden-headed lion tamarins (Leontopithecus 
chrysomelas) in Brazil; none of the animals had skin lesions. 
Similarly, none of the 74 bats from which isolates of M. ves- 
pertilionis were recovered in this study showed visible signs 
of dermatitis at the time they were sampled. Thus, there is no 
current evidence that M. vespertilionis acts as a pathogen; 
instead, the yeast is likely a normal component of the skin 
mycobiome of bats. 


In this study, M. vespertilionis was isolated from nine species of 
bats in seven U.S. states. With the exception of Lasionycteris 
noctivagans, all of these bats were species of Myotis, a diverse 
and widely distributed genus in North America. This, coupled 
with the detection of the fungus in both the eastern and western 
United States (Fig. 2), may indicate that the yeast is more wide- 
spread than our survey indicates. Myotis also has representa- 
tives throughout much of South America, Eurasia, Africa, and 
Australia, and additional sampling is necessary to determine 
whether M. vespertilionis may have a more global distribution. 
If this is shown to be the case, M. vespertilionis could prove to 
be an important species in which to study Malassezia genetics 


and evolution. All other known species of Malassezia occur 
on humans and domestic animals, which have been trans- 
ported across the world. This, combined with the possibility of 
interspecific recombination (Cabafes et al. 2007, Castella et 
al. 2014) when different species come into contact with one 
another, makes it difficultto understand historic patterns of geo- 
graphic distribution, variation, and genetic exchange between 
strains. Bat species, however, have not been subjected to such 
extensive and recent human-assisted global movements, nor 
do they frequently come into close contact with other animal 
hosts that might facilitate genetic transfer between different 
species of Malassezia. Thus, if M. vespertilionis is found on 
other continents, genetic distinctions between strains may still 
be intact and provide important opportunities for future research 
into the evolution of these fungi and the role that humans have 
played in shaping the genetic structure and pathogenicity of 
Malassezia species. 
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